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ABSTRACT 

The measures of mechanical ahgnment were obtained for both prolate and 
oblate grains when their temperature is comparable with grain kinetic energy 
^ ' devided by k, the Boltzmann constant. For such grains, the alignment of 

angular momentum, J, with the axis of maximal inertia, a, is only partial. This 



0^ 






substantially alters the alignment as compared with the results in Lazarian 
(1995) and Roberge, Hanany & Messinger (1996) obtained on the assumption 
>• ■ of perfect alignment. We also describe the Gold alignment when the Barnett 

dissipation is suppressed and derive an analytical expression which relates the 
measure of alignment with parameters of grain nonsphericity and the direction 
Q ■ of the gas - grain drift. This solution provides the lower limit for the alignment 









measure, while the upper limit is given by the analytics derived in Lazarian 
(1994). Using results of a recent study of incomplete internal relaxation in 
CL|' Lazarian & Roberge (1996), we find measures of alignment for the whole 

O ' range of ratios of grain rotational energy to k over Tg, where Tg is the grain 

-*-» ' temperature. To describe alignment for mildly supersonic drifts, we suggest 

c^ ■ an analytical approach which provides good correspondence with the results of 

direct numerical simulations in Roberge, Hanany & Messinger (1995). We also 
^ ■ extend our approach to account for the simultaneous action of the Gold and 

d • Davis- Greenstein mechanisms. 



Subject headings: dust, extinction — ISM, clouds — ISM, polarization 



1. Introduction 

Understanding the observed alignment of the ISM grains is yet unsolved astrophysical 
problem (see Hildebrand 1988, Whittet 1992, Goodman et al. 1995, Roberge 1996). This 
limits the use of polarimetry data for studying interstellar magnetic fields. 
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The mechanism of mechanical ahgnment of thermally rotating grains was pioneered 
by Gold (1951) the same year that the classical paper introducing the paramagnetic 
alignment by Davis & Greenstein (1951) was published. Originally, Gold suggested that 
grain alignment arises from cloud - cloud collisions (Gold 1951, 1952), but it was shown 
in Davis (1955), that such collisions can align only an insignificant fraction of interstellar 
grains. Therefore a further study of the mechanism was devoted mainly to the alignment 
in the vicinity of bright sources, where radiation pressure can drive grains to supersonic 
velocities (see Purcell 1969, Aitken et al. 1995). A new important idea, put forward by 
Roberge & Hanany (1990), that grains can be aligned by ambipolar diffusion made the 
Gold alignment more promising (see also Roberge, Hanany & Messinger 1995 hereafter 
RHM). Our study in Lazarian (1994, henceforth Paper I, and 1995a) showed that the role 
of mechanical processes had been underestimated. It was found that pervasive MHD wavesQ 
can produce grain alignment even within ideal MHD. 

The goal of this paper is to provide an analytical quantitative description of the 
mechanical alignment of thermally rotating grains, refered to as "Gold alignment" . We 
distinguish between the Gold and mechanical alignment of suprathermally rotating grains 
(see Lazarian 1995a,b, Lazarian & Efroimsky 1996, Lazarian, Efroimsky & Ozik 1996). 

Although the expressions for the measure of alignment corresponding to Gold alignment 
were obtained in Paper I, the shortcoming of this study was that only perfect alignment 
of angular momentum, J, with the the principal axis of the maximal moment of inertia, 
a, (henceforth the axis of major inertia) was considered. Although this regime is valid for 
highly supersonic motions discussed in Lazarian (Paper I), the wealth of the ISM conditions 
presents us with a wide range of other options. 

It is well known from theoretical mechanics, that internal dissipation of energy cannot 
change grain angular momentum and, for a grain with fixed angular momentum, the 
minimal energy corresponds to rotation about the axis of major inertia. For suprathermally 
rotating grains, efficient internal dissipation of energy leading to nearly perfect alignment 
of angular momentum with the grain axis of major inertia was discovered by Purcell (1979) 
(see also Spitzer & McGlynn 1979). For thermally rotating grains, the alignment is only 
partial (Paper I, Lazarian & Roberge 1997, henceforth LR97). 

To relate polarimetry observations, that are influenced by the alignment of the long 
grain axis, and the theory, that deals with the alignment of angular momentum, J, one 



^Here and further on in the paper when speaking about MHD waves we mean both Alfvenic and 
magnetosonic waves. 
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has to describe the alignment of J not only in respect to magnetic fieldQ, but also in 
respect to grain axes. These can be called external and internal alignment, respectively. 
Internal alignment arises both from the difference in grain moments of inertia ( "Maxwellian 
alignment") and due to the Barnett relaxation^ ("Barnett alignment"). 

Below we use the statistical approach introduced in LR97. The gist of it is to describe 
deviations of J from the axis of major inertia, a, as thermal fluctuations from the position 
of static equilibrium when the time of internal relaxation is much shorter than that of the 
gaseous damping. By introducing the grain rotational temperature, which depends on both 
the temperature of the ambient gas and on the drift velocity, we show that, when the ratio 
of this temperature to the temperature of grain material tends to infinity, J is coupled with 
a, and this justifies our approach in Paper I. For finite ratios however, the deviations of J 
become important for the alignment. 

To describe Gold alignment for mildly supersonic velocities, we have to account for 
the rms motions of gas atoms. To do this, we assume that grains are subjected to a 
superposition of fluxes within a small angle to the direction of the initial flux and the value 
of this angle depends on the ratio of the flux velocity to the mean velocity of thermal 
motions. This enables us to obtain an asymptotic solution which provides fair agreement 
with the numerical calculations in RHM for the case of perfect coupling between J and 
a. Within the same model of perfect coupling the agreement with numerics in RHM is 
obtained for the joint action of Gold and Davis-Greenstein mechanisms. A comparison 
between numerics and analytics for the general case of incomplete internal alignment is 
planned, and we anxiously follow the progress of numerics in this direction. 

The structure of the paper is as follows. First of all, we remind our reader of the 
concept of internal relaxation and describe how grain drift influences internal alignment 
(Sect. 2). Then in Sect. 3, we obtain analytical solutions for the alignment measure when 
internal relaxation is negligible and compare these results with the other extreme case, 
namely, when J is perfectly coupled with a. The treatment of the problem for an arbitrary 
degree of internal alignment is given in Sect. 4. In Sect. 5, we compare our results with 
numerics and argue that our analytical approach provides an adequate description of the 
alignment both for mildly supersonic velocities and for Gold and Davis-Greenstein processes 



•^Magnetic field acts as the axis of alignment because the time-scale of Larmor precession is usually much 
shorter than the time-scale of alignment. We remind the reader that the rapid Larmor precession arises from 
magnetic moments of grains; those moments are caused by the Barnett effect (Dolginov & Mytrophanov 
1976). 

''As a rule, the Barnett relaxation dominates internal relaxation (see Purcell 1979). 
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acting together. A short summary of results is presented in Sect. 6. 



2. Internal alignment 

As pointed out above, to relate theory and observations it is essential to describe the 
alignment of grain long axes in respect to magnetic fieldQ. This alignment is determined 
uniquely by the alignment of grain angular momentum only if the angular momentum is 
coupled with the axis of major inertia, a. Further we study a general problem when the 
internal alignment is partial. 

Theoretical studies deal with the measure of alignment of angular momentum J|^ 

Q, = ^{cos'9^-^) , (1) 

where 6i is the angle between J and the direction of magnetic field, while polarimetry 
provides us with the data on the Rayleigh reduction factor (Greenberg 1968) 

R=l{cosH)-l , (2) 

where 6 is the angle between the direction of magnetic field and the symmetry axis of a 
spheroid approximating the grain. The question how to relate these two quantities was on 
the astrophysical agenda from the very beginning of research in the area. First in Jones 
& Spitzer (1967), it was assumed that the distribution function of J in the grain reference 
frame is independent of the alignment of J in respect to magnetic field. This was a natural 
assumption to start with. Later, Spitzer (1978) showed that additional ahgnment of J in 
respect to the axis of major inertia should be present due to paramagnetic relaxation in the 
external magnetic field. This effect can be called Spitzer relaxation to distinguish it from 
the Barnett relaxation discovered by Purcell (1979). 

Here we assume that the ratio of the gas damping time to the time of internal 
relaxation is much greater than unity (see estimates in Roberge, DeGraff & Flaherty 



^In a special case, e.g. in the atmospheres of comets, the ahgnment can happen on the time-scale less 
than that of precession. In such situations, the alignment should be defined in respect to the flow, which is 
similar to the alignment when the flow is directed along magnetic field and therefore we do not discuss it 
separately. 

^Our notations are different from the ones used in Paper I. In an attempt to insure that the system of 
notations used by different authors is universal we adopt the system of notations suggested in review by 
Roberge (1996). 
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1993). In this case, the deviations of J from the axis of major inertia, a can be described 
thermodynamically (LR97). 

For a grain with fixed angular momentum J = | J| and with components of the moment 
of inertia related as I^ < ly < h, the kinetic energy is 

J2 



E = — 
2 



;i/4 - l//y) sin" e sin" d + (V^y - l/h) sin" 6 + 1//, , (3) 



where ^ is the azimuthal angle in the x — y plane, and 6 is the angle between J and 2;-axis. 
Although the equilibrium position of the grain with fixed J^ corresponds to ^ = 0, thermal 
fluctuations cause deviations from it. These deviations can be described by the Boltzmann 
factor exp{-E/{kT,)} (LR97). 

A complex precession of J due to these fluctuations modifies both the grain interaction 
with the gaseous flow and the dichroic absorption. However, it seems reasonable to assume 
that the behaviour of such an asymmetric grain can be approximated by the behaviour of 
a spheroidal grain with the moment of inertia somewhere between I^ and ly (further on 
we denote this mean moment by /_l). Here we also adopt this approximation and hope to 
compare the results for spheroidal and irregular grains elsewhere. 

For a spheroidal grain with termprature Tg, Eq. (^ can be simplified, and the 
distribution function of the angular momentum in the grain reference frame is (LR97) 



/te(6'2) = const X sin6'2exp 



^ {ih-l)sm'92 + l} 



2LkT. 



(4) 



where 62 is the angle between J and the rotational symmetry axis of the spheroid, h = Iz/I± 
and the value of J is different for different grains of the ensemble. The fact that the 
probability of a particular angular momentum depends on angle 62 complicates the study. 
However, computations in LR97 show that with a sufficient degree of accuracy it is possible 
to substitute the rms value of J for a Maxwellian angular momentum distribution in Eq. (1) 
to get an approximate measure of alignment for an ensemble of grains with given rotational 
temperature. Therefore as the zero approximation we use 

(5) 









J2 = 


kT,sL{2/h + l) 


(Landau 


& Lifshitz 


1980), 


where 


T„ + Ts 1 mu^ 

— ^ H 

2 3 2k 



(6) 

m and u are the mass and velocity of bombarding atoms, respectively, and Tg is the gas 
temperature. We expect that numerical simulations by Roberge will determine the accuracy 
of the approximation adopted. 
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For hypersonic drifts in diffuse medium discussed in Paper I, T^s ^ T^ and the 
ahgnment of J in the grain reference frame is close to being perfect. Indeed, if the drift 
velocity of grains exceeds the velocity of atoms in diffuse clouds 3 times, Ts/T^fi ~ 0.1. The 
corresponding distribution of J (see Eq. (^) has a peak near the axis of major inertia (see 
Fig. 1) and can be approximated by a delta- function. This, however, may not be true for a 
different grain environment. For instance in molecular clouds, gas and grains have similar 
temperatures and, for mildly supersonic drift velocities, T^s is of the same order as T^. 

In short, perfect alignment between J and a does not seem to be universally applicable 
to interstellar grains. Below we describe the Gold alignment when this constraint is lifted. 



3. Gold alignment for J not parallel to a 

3.1. Analytics for the Rayleigh reduction factor in the absence of internal 

relaxation 

The distribution of angular momentum can be characterized by function f{n,J), 
where n is the number of grain-atomic collisions. In general, the direction of J should be 
defined by angles 6i and ipi in the "gas reference frame" and by 62 and (f2 in the "grain 
reference frame" (see Fig. 2). Angle ipi describes the precession of J about magnetic field 
and angle ip2 describes the precession of the spheroid's axis of rotational symmetry about 
J. Henceforth grains will be approximated by spheroids with semiaxes a and b. 

As the change of grain angular momentum in the course of an individual collision is 

small, the alignment of J can be described by the Fokker-Planck equation (see Reichl 1980, 

Roberge, DeGraff & Flaherty 1993). Then following Dolginov & Mytrophanov (1976), we 

can write 

■9/(x,n) 2 g/(x,n) 2 gV(x,n) 

where x is a vector in the phase space with coordinates J, cos 6*1, cos 6*2, fi and ip2, and the 
coefficients Oj and bik obtained on the assumption of supersonic grain drift can be found 
in Appendix A. The solution of Eq. (^ in the limit of hypersonic drift is (Dolginov & 
Mytrophanov 1976) 

constg f J'^{l + gcos'^92 + scos'^ei)\ ^^^ 

J{j,cos0i,cos02,n) = — :——exp' 



n^/2 \ 2nb'^p'^{l + g + s 

where 



2 (p2) - (p2) ^y^ 
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is the external flux anisotropy and 

' = -liT ^''^ 

is the grain non-sphericity. Note, that (p^) and {pD are the averaged squared momentum 
and its Zi component (see Fig. 2) transferred to a grain in an individual collision. Both g 
and s can vary from —0.5 to oo. It is easy to see that g = —0.5 corresponds to flakes and 
(jf — > oo to needles; s = —0.5 corresponds to a flux perpendicular to magnetic field, while 
s ^ oo to a flux parallel to the field. The fluxes are measured in the grain reference frame, 
and therefore a gaseous flux with the velocity u is equivalent to grain drift with the velocity 
—u in respect to the ambient gas. 

Spherical grains {g = 0) correspond to a = 6, while isotropic fluxes {s = 0) to 
(p2) = 3 {pD- We expect changes in grain alignment when s = and/or g = 0. To have 
a picture which is easy to visualize, we refer to fluxes (drifts) corresponding to s < as 
"fluxes (drifts) at large angles to magnetic field" , and to those corresponding to s > as 
"fluxes (drifts) at small angles to magnetic field". Note, that "small" angles lie within 
the interval [0, arccos(l/-\/3)], while "large" angles in [arccos(l/A/3), 7r/2], if we limit our 
discussion to the first quadrant. Obviously, oblate and prolate grains correspond to g < 
and g > 0, respectively. 

Calculations in Dolginov & Mytrophanov (1976) show that the solution given by 
Eq. (1) is accurate up to 0.251(751 order terms for |s| < 1 and l^^l < 1. The accuracy becomes 
of the order of s"^ when s —>■ oo and \g\ < 0, whereas it is g~'^ when g —>■ oo and s < 0. If 
both g and s are large, the accuracy is of the order of max[(yf~^, s~^]. 

Angle 6 can be found from simple geometric considerations (see Fig. 2) 

cos6' = cos6'i cos6'2 + sin6'i sin6'2 cos(v9i — (/)2) • (H) 

Averaging out the dependencies on ipi and ip2 provides (see Davis & Greenstein 1951 
eq. 108): 

(cos2^)^ = 0.5(l-cos2^i-cos2e2 + 3cos2eicos2e2) ■ (12) 

To calculate the Rayleigh reduction factor given by Eq. (0), it is necessary to find {cos'^9). 
After averaging Eq. (^ over J the required distribution function can be found to be 

WG{eue2) = C{g,s){l + gcos^e2 + scos^ei)-l , (13) 

where C{g, s) is the normalization constant for any fixed s and g. Then 

2 m _ II cos^d Wg singising2dgidg2 . ^ 

^^°^ ^~ JJWg sin 01 sin 02(161(162 ' ^ ' 



The necessary calculations are performed in Appendix B, and here we present the final 
expression 



where 



{cos^e) = — {l + g + s + gs-C{g,s)^l+g + s) , (15) 



C{g,s) 



y/9S Urctan ^j^j , ^s > 



(16) 



^s| (arctanhW J^j , ^s < 



which is valid for all possible values of s and g. This expression can be compared with the 
analytics found for the case of strong relaxation (see Appendix C). 

The corresponding Rayleigh reduction factor is 



This measure enters the formulae for intensity of polarized radiation due to dichroic 
absorption. Equation ([T7| ) encompasses a variety of circumstances, which are explored 
below. 



3.2. Comparison with Paper I 

Here we discuss the alignment for various values of s and g. When internal dissipation 
aligns J and the axis of major inertia perfectly, Eqs. (|12|), (Q) and (|I]) give the following 
relation between the Rayleigh reduction factor and the measure of alignment of angular 
momentum, aj, used in Paper I 

^ ^ f o-j, for oblate grains , 
I — O.ScTj, for prolate grains . 

It is easy to see that prolate and oblate grains produce polarization of the same sign despite 
the fact, that the Rayleigh reduction factor has opposite signs in these two cases. This 
becomes clear if one recalls that oblate grains aligne their short axes in respect to the 
magnetic field, while prolate grains align their long axes. 

First, consider oblate grains {g < 0) subjected to drift at large angles to magnetic field 
(s < 0). The corresponding measure of alignment is shown in Fig. 3. A comparison with 
fig. 3 in Paper I reveals the decrease of alignment associated with suppression of internal 
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dissipation. For instance, for the drift perpendicular to magnetic field lines, this measure for 
flake-like grains is 0.25 if the Barnett alignment is absent, while it reaches unity (complete 
alignment) if the internal dissipation is efficient (see Fig. 4). 

The alignment measure for oblate grains drifting at small angles to magnetic field 
[s > 0) is shown in Fig. 6. The comparison with fig. 4 in Paper I testifies an order of 
magnitude decrease of the measure. A cross-section of the plot for g = —0.5 (see Fig. 5.) 
shows that fiakes are only marginally aligned. This is in contrast with the case of intense 
internal dissipation, when fiakes are well aligned and the measure of alignment approaches 
—0.5, if the fiux and magnetic field directions coincide (Lazarian 1994a). Due to the 
symmetry of the distribution function given by Eq. (^ under a simultaneous interchange 
cos^i ^^ cos ^2 and s ^^ g, we can claim that needles are marginally aligned if grain drift 
is perpendicular to the field, i.e. u±B. Indeed, the corresponding measure of alignment 
for prolate grains is negligible (see Fig. 7). In breif, if the internal dissipation is negligible, 
prolate and oblate grains are marginally aligned for drifts at angles close to 7r/2 and 0, 
respectively. 

When both s and g are positive, the corresponding measure of alignment is shown in 
Fig. 8. One can see that this measure tends to +0.25 when s or/and g tend to infinity. 
In terms of aj (see Eq. ([T8|)) this is equivalent to crj = —0.5 obtained in Paper I. This 
correspondence has a simple explanation. Indeed, for a needle the angular momentum 
should be directed along the axis of major inertia even in the absence of internal dissipation. 

The requirement of large s and g places stringent constrains on the efficiency 
of alignment for typical ISM conditions. Indeed, let grain drift velocity components 
perpendicular and parallel to magnetic field be u± and uy, respectively, then if u± is much 
greater than the rms velocity of gas atoms Vrms 

where w ~ u±/u\\. If u± -C Vrms, than w ~ Vrms/v\\ should be used in Eq. (p!9|). In any 
case, w is unlikely to be less than 0.1, and therefore s is not likely to be greater than 100. 
We neither believe that the axis ratio of a typical prolate grain exceeds 10, therefore g is 
likely to be less than 50. The measure of alignment for w G [0.2,0.7] and the axis ratio 
y = a/b G [1, 10] is shown in Fig. 9. The cross-sections of the plot for y = 5 and y = 10 are 
shown in Fig. 10. 

To summarize, in the absence of internal dissipation, the measure of alignment drops. 
Oblate grains are most aligned for u_LB, and prolate grains are most aligned for u||B. 
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4. Generalized problem 

Up to now two extreme cases were discussed: strong internal dissipation, when the 
angular momentum is coupled to the axis of major inertia, and weak internal dissipation, 
when the residual alignment of angular momentum is due to the differences between the 
maximal and minimal rotational inertia. In both cases analytical solutions were found. 
These solutions provide the upper and lower bounds for the measure of alignment. For 
instance. Fig. 11 shows these two bounds for u±B when the alignment is caused by 
ambipolar diffusion. A conspicuous feature of this particular figure is that the two plots are 
different for ideal spheres, i.e. when no axis alignment is expected. Equation (^ testifies 
that for /i — i> 1, all positions of J become equally probable. This fact was ignored in our 
simplified approach adopted in Paper 1. Fig. 11 also testifies, that the internal dissipation 
strongly influences grain alignment, as the spread in Rayleigh reduction factors for the two 
extreme cases is wide. 

To account for the incomplete alignment of J in the grain reference frame, one has to 
incorporate internal dissipation in the Fokker-Planck equation. An analytical study in this 
case seems formidable, and a numerical approach, e.g. similar to the one used by Roberge, 
DeGraff & Floherty (1993) may be advantageous. For such a study both the analytical 
solutions obtained above and those derive in Paper 1 should serve as benchmarks. 

A less rigorous, but less laborious way to account for the incomplete alignment is 
to follow Jones & Spitzer (1967). Let the distribution function of J in the presence of 
internal dissipation, VFgd, be the product of the distribution functions Wq and Wy) given 
by Eqs (|1^) and (^), respectively. Then after expressing h in terms of g and J^ in terms of 
Teff using Eqs (|^) and (||), respectively, we get 

l^GD - const X sin 02 (l + g cos^ 62 + 3 cos^ ^i)' 2 exp '^^ ^ ^ -\ . (20) 

\ -'s 9 ' ^ ^ ) 

Using H^GD in Eq. ([I^ ) we calculated (cos^ ff) and the Rayleigh reduction factor for 
Tcs/Ts equal to 100, 10 and 1.1. For example. Fig. 12 shows the Rayleigh reduction factor 
as a function of grain oblateness. The comparison between Figs 12 and 11 shows that while 
for efficient internal relaxation corresponding to T^q/T^ > 100 the approximation of perfect 
coupling is appropriate, in the case of Trot ~ Ts, this is no longer true, and the analytics 
disregarding internal relaxation provide a better fit. Similar conclusions are valid for other 
values of s and g with the exception oi g —>■ 0. In this case, as discussed earlier, there is no 
coupling between the angular momentum and the axis of greatest inertia. 

The alignment of prolate grains by ambipolar diffusion is marginal and the comparison 
between Figs 13 and 7 shows that the internal dissipation does not change it much. In 
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contrast, Figs 14 and 5 show that the internal dissipation drastically changes the alignment 
of oblate grains due to the radiation pressure. If this alignment is marginal when the 
internal relaxation is suppressed (see Fig. 5), it becomes substantial as soon as the internal 
relaxation is present. This difference in the susceptibility to internal dissipation is easy to 
understand if one recalls that without internal dissipation J is only marginally aligned with 
the axis of major inertia for oblate grains and the alignment is substantial in the case of 
prolate grains. Therefore Fig. 14 shows a marginal difference in grain alignment for different 
values of internal dissipation as grains become sufficiently prolate. 

To summarize, our results show that Gold alignment is modified by internal relaxation, 
and to predict the alignment accurately it is essential to estimate the ratio of the effective 
rotational grain temperature to grain material temperature. Small temperature ratios 
usually occur when the grain drift and the thermal velocities are comparable.^ Then it is 
essential to account for the fact, that, in the grain reference frame, atoms move at the drift 
velocity, which is modified by thermal motion. In other words, grain - gas collisions formally 
correspond to a range of values of s. To account for this effect, we suggest to integrate 
over the corresponding range of s (see Lazarian 1995a). Since there is no direct numerical 
calculations of the measure of alignment for the case of incomplete internal relaxation, in 
what follows we compare our results only with numerics for perfect coupling between J and 
the axis of major inertia. 



5. Comparison with RHM 

A comprehensive numerical study of the Gold alignment for perfect Barnett relaxation 
was done in RHM. The authors presented a detailed study of grain alignment for a range 
of drift velocities starting from subsonic ones. Paramagnetic relaxation was also included 
in their model. Here we briefiy discuss how to improve our model to include both subsonic 
drift and paramagnetic relaxation. 



5.1. Subsonic drift 

The model adopted here assumes that the grain drift is essentially hypersonic. In 
other words, we have ignored the rms velocity of gaseous atoms as compared with the drift 



^We disregard a rather artificial case of hot grains drifting in cold gas. 



-12- 



velocity and assumed that atoms hit the grain from one direction defined by 



arcsm 



"^ + "^ (21) 



\ul + u'^ + ul 



where Ui, i= x, y, z, are the components of atom drift velocity in the grain reference frame. 
Obviously, the anisotropy parameter s can be expressed as a function of (p 

s = ^-^'^f^ (22) 

1 - C0S2 ^ ' 

and hence the Rayleigh reduction too. 

When the drift is subsonic, in the grain reference frame atoms are viewed as approaching 
from various directions, any particular atom at angle 



arcsm , 



(Wx + ^x)^+(Wy + fy)^ f^^. 

where Vi, i= x, y, z, are the components of rms velocity of the atom. According to Eq. (^31) , 
angle varies from atom to atom due to variations in Vi. To obtain the statistics of R we 
subdivide atoms into subgroups of atoms having the same f i.0 For each of the subgroups, 
is constant, hence the measure of alignment can be obtained by averaging R{ui,Vi) over a 
Maxwell-Boltzmann distribution of fj. 

Elsewhere we hope to test the applicability of such an approach for a considerable 
range of drift velocities by comparing our predictions with direct numerical simulations. At 
the moment, we want to show that our estimates are in a reasonable agreement with the 
data presented in RHM. 

RHM assume that the angular momentum is perfectly coupled with the axis of major 
inertia, which corresponds to Ts/Tcs = in our model. The effect of the spread of atom 
velocities in the grain reference frame should not depend upon the position of the angular 
momentum in respect to the axis of the major inertia. Indeed, the velocities of particles 
on the grain surface are much smaller than the velocities of striking atoms and therefore 
if we obtain the correspondence between our and RHM predictions we may hope that our 
treatment is applicable when the internal alignment is incomplete. 

Figure 8 in RHM shows the Rayleigh reduction factor for oblate grains drifting with 
different velocities in respect to gas. The calculations are made up to Mach number 4, but 



^As the time of alignment is much shorter than the damping time (see Paper I), each of the subgroups 
may by itself cause alignment. 
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the saturation of the alignment is obvious from this plot. Therefore we will consider that 
these values correspond to hypersonic velocities and define -R(0o); where 0o in the case of 
ambipolar diffusion studied in RHM is 7r/2. For our simplified estimates we observe that 
for large Mach numbers M , Eq. (PSf) provides 



arcsin (1 — 0.5 M 



2 2 24 

where the first term corresponds to (pQ and the rest can be interpreted as 

1 



{Rh 



rfpo+Sfp dR ds 



where 



R{(J)q) = -3.5 -3^ + 3 Jl + 2^(1+^) arcsin 



1 



2v/TT^ 



Therefore 



(24) 

(25) 
(26) 



dR 

d^ 



3 + 3(7 



s=-0.5 



v/T+2^v/3 + 4^ 



+ a/3 + Ag {A + Q g) arcsin 




(27) 



and 



(2^ 



ds 6 cos a; sin a; (2 — 6 cos^ x) cos x sin x 
dx l-cos^a; (l-cos^x)^ 

with ds/dx ~ 0.032 when x = (p = 7t/2. 

For instance, RHM find that at large Mach numbers the Rayleigh reduction factor 
for grains with the axis ratio 0.5 is ~ 0.25. This axis ratio corresponds to g = 1.7, when 
pluged into Eq. (^) it gives {R)(f, ~ 0.21 for M = 2. This is comparable with the value 
^ 0.19 that follows from fig. 8 in RHM. Similarly for the axis ratio 0.25 corresponding 
to g = 7.5, {R)(f, ~ 0.27, for M = 2, which is of the same order, as the result in RHM 
{^ O.26).0 As the values of the Rayleigh reduction factor obtained in RHM for high Mach 
numbers essentially correspond to the values obtained in the analytical treatment in Paper I 
(see fig. 12 in RHM) it is possible to see, that a purely analytical treatment is appropriate 
at least for some values of subsonic drift velocities. The entire range of velocities will be 
treated elsewhere. 



^Again we take the values of the Rayleigh reduction factor obtained for M 
in our formulae. For the axis ratio 0.25 this value « 0.39. 



4 in RHM and substitute 
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Our estimates above were obtained using the analytical solutions obtained in Paper I 
for perfect coupling of J with the axis of major inertia. Evidently our approach is 
also applicable to describe alignment at low Mach numbers using the analytics and 
"semi-analytics" obtained in Sections 3 and 4 respectively. We are looking forward to the 
progress in numerical techniques to be able to compare our predictions with the results of 
direct numerical simulations. 



5.2. Mechanical and paramagnetic alignment 

So far we have talked only about mechanical alignment and completely disregarded the 
paramagnetic one. This is justifiable only if the Davis- Greenstein alignment is negligible. 
In general, the Davis- Greenstein alignment must be accounted for. To do this we propose a 
simple formular. 

If we denote the Rayleigh reduction factor of the Gold alignment by Rq, the measure of 
the internal alignment by Qx, and that of Davis-Greenstein process by Rdg^ it is possible 
to estimate the measure of the overall alignment as 

Qxcr g + RgRdg + Qx RpG 
Qx + '^RgRdg 

To obtain the expression above we used the expression for J alignment given by eq. 45 in 
Lazarian (1995a) and the approximation 



^ ~ ^^ 7^2 — I on J-, • \'^^l 



Ri ^ Qjii) X Qx , (30) 

where Ri is the Rayleigh reduction factor obtained for the i mechanism acting alone, while 
Qj[i) is the measure of J alignment relatively to magnetic field. The latter approximation 
follows from spherical trigonometry (see Eq. (|T^)). Indeed, if two processes are independent, 
then 

(cos^ 01 cos^ 62) ^ (cos^ ^1) (cos^ 62) . (31) 

The approximation above was proved to be sufficiently accurate for the Davis-Greenstein 
process in Lazarian (1995a), but may be much less accurate for the Gold alignment. In any 
case, we treat Eq. (^) only as a conjecture to be tested in future. 

A study of the simultaneous action of the paramagnetic and mechanical alignment has 
been done recently in RHM for T^/Tg = 0. In this case, 0"^ = 1 and Eq. (^) reduces to the 
one derived in Lazarian (1995a). 

First of all, to provide the comparison we must define the measure of alignment for 
the Davis-Greenstein process when grains are subjected to a supersonic fiow. Naturally, as 
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grains drift faster the rate at which atoms arrive to their surface increases, and therefore 
the gaseous damping time decreases. For our rough estimate of the Davis-Greenstein 
ahgnment subjected to the supersonic flow, we simply substitute overall velocity instead of 
the thermal one into the expressions for the diffusion coefficients for gaseous bombardment 
(see Roberge, DeGraff & Flaherty 1993). This changes both Tcs and the ratio of the gaseous 
and magnetic damping times 6i (see eq. 25 in Lazarian 1995). The ffrst change is not 
relevant here because to compare our estimates with the calculations in RHM, we assume 
Ts/Tcs = 0. On the contrary, the second change rescales 6i calculated for the non-drifting 
grain (this value is shown in figs 4 and 9 in RHM). Let the rescaled value be 6/M. To 
obtain ctg, we use the zero approximation in Lazarian (1995c). 

For spherical grains when 6i = 10, RHM obtain a ^ 0.40 (see fig. 4 in RHM). Our 
approximation gives a ^ 0.38. It is easy to see from fig. 9 in RHM that for 6i = 1, axis ratio 
0.5 gives a ^ 0.32, whereas axis ratio 0.25 gives a ^ 0.48. Our approach provides a ^ 0.29 
and 0.46, respectively. This approximate correspondence let us hope that our simplified 
analytical treatment reproduces essential features of the alignment. 



6. Discussion 

In short, we have shown that the alignment of J in the grain reference frame, that arises 
from the difference between the grain material and rotational temperatures, is essential for 
the Gold alignment. When velocities of grain drift are hypersonic, we can assume perfect 
coupling of J and the axis of major inertia. However, this assumption fails for mildly 
supersonic drifts when the gas and grain temperatures are comparable. Such conditions 
are expected, e.g. in molecular clouds undergoing ambipolar diffusion. Then incomplete 
alignment between J and the axis of major inertia should be accounted for. 

Our results also show that the effect of incomplete relaxation is more vivid for oblate 
grains than for prolate ones. This is a consequence of the fact, that for sufficiently prolate 
grains the alignment of J in respect to grain axis of major inertia is manifest even without 
internal relaxation. 

The analytical results obtained above and those derived in Paper I provide the lower 
and upper bounds for Gold alignment. In a general case of incomplete internal alignment, 
we suggested a semi-analytical approach. 

To provide a quantitative description of Gold alignment for drift velocities comparable 
with the thermal velocities of gaseous atoms, we suggested an expansion of Rayleigh 
reduction factor in a series over the drift Mach number and obtained a fair correspondence 
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with the results of direct numerical computations. 
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A. CoefRcients of the Fokker-PIanck equation 

The coefficients a, and 6,^ of the Fokker-PIanck equation are determined by the change 
of the grain angular momentum due to grain - atom collision 

/l. SAxi 1. (9Aa;i .\ ,.^, 

/l d AXm I . dAxm. 

A xi — h - A x„ 



5ifc 



\2 dxi 2 dxm 

+ I A ¥^™|^ - Ax„\ , (A2) 

2 O Aifm I 

{Axi A Xk) , (A3) 



where m = 1,2, k,i = 0, 1, 2, and the angular brackets denote averaging over atom impacts 
over the grain surface. The quantities Axj and Acpi and the corresponding a^ and bik were 
calculated in Dolginov & Mytrophanov (1976) using the following equations: 

A Xj = A{ej ■ J) - (e^ ■ J) A J , 

Axo = j ■ AJ , 

Ay,, = ((e,xj).AJ(J(l-a;2))-i , 

where j = 1,2, ei = H/|H| is a unit vector along the magnetic field, e2 = a/|a| is a unit 
vector along Zi-axis of the grain, j = J/| J| is a unit vector along J, and AJ = r x p with r 
and p used for the vector to the point of atomic impact and atom momentum, respectively. 

To find coefficients a^ {k = 0, 1, 2) and bik {h A; = 0, 1, 2) (see Dolginov & Mitrophanov 
1976, eq. 5) one has to substitute Axj, Axq, Aipj into Eqs (|A1|) , (^), and (|A3|) and 
perform the necessary averaging. Apart from averaging over the surface area exposed to 
the fiux, one has to average over the angles of precession of of J about e2 and m. 



-17- 



B. Computation of integrals 

Equation (jM]) involves integration of Eq. (|13|) with the distribution function Wg{Oi, O2). 
As a result, we get 

(cos^ 6) = 0.5 (1 - (cos^ Oi) - (cos^ 62) + 3(cos^ Oi cos^ ^2)) , (Bl) 

where the explicit expression for the first term is 

/ 2/)\ r-f \ r'^ Aa r'^ cos^gisin^isingg ^^ ,^^. 

{cos 61) = C{g,s) dOi -— — — 2/3^3/2 ^^2 • B2 

Jo JQ [1 + scos'^Oi+ gcos^ O2Y' 

After integrating over 62 and the change of variables, sin^i = x, we have 

(cos^ e,) = C{g, s) / , (B3) 

where 



Using the identity 
Equation (p^) reduces to 



-(l + sx^)-- = x2 , (B5) 



(cos^gi)=C(^,s)-/ -- , (B6) 

s JO V 1 + sx^ + g s 

which can be calculated [Gradshtein & Ryzhik 1965, 2.271(4)] to give 

r 1 In Vs±^±s ^>o 

^1 = / , = (B7) 

•^° V(l+^) + ^^' [^arcsin^^, s<0 . 

For uniform representation, one can also use inverse hyperbolic function for s < 0, namely 



1 / s 
ii = —= arcsinhJ , (BS) 

V^ yi + g 



where 



arcsinh z = ln(z + y/z'^ + 1) = -arcsin (iz). (B9) 



I 



To find C(g, s) = 12 ^, we calculate 



'' - L [l + sx^)Jl + sx^ + g ■ ^^^°^ 
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By substituting u = x'^ + s into Eq. (|B10|) and evaluating the resulting integral (Gradshtein 
& Ryzhik 1965, 2.224(5)), we have 



t2 = { (Bll) 



^ ^ In v+^+9+v-3^ < 



_|_ arctany^3r, ^^^ > q 



Using the inverse hyperbolic function 



1 1 + z 
arctanh 2; = - In = arctan [iz] , (B12) 



the expression for 12 can be rewritten as 



1 / —gs 

^2 = ^^^= arctanhw , sg <Q . (B13) 

V=^ V 1+^+^ 



Finally, we obtain, 
for s < and (7 < 0, 



for s < 0, and (? > 0, 



for s > 0, and (7 < 0, 



for s > 0, and (7 > 0, 



\/^ arcsin./— 77- 1 
s arctan./ -. ,^ , s 



,/g arcsin./— 77- 1 

(cos2^i) = ^^^^ ^^±L-i , (B15) 

s arctanh,/— -. ,^^, s 



-q arcsinh./TT- 1 

(cos^^,)= ' ^ I ^]V -- , (B16) 

s arctanh,/— -. ,^ , s 



v^ arcsinh^^ 1 

(cos ^1) = ^^= - - , (B17) 

s arctan^^^ s 

which covers all cases. Similarly, for the second term in Eq. (pl|) we obtain 



^^3,s)^j^ arcsiny^-Y^, ^ < , 
(cos' 62) = { (B18) 

^^^' ^) V5 arcsinh^^, ^ > . 

The last term in Eq. ( pi|) is 

/Q 2 /J 2n\ o^. i-n/2 ^n/2 cpg^ g^ cos' ^2 sin ^1 siu 62 d^i d^a .^^„. 

3cos ^icos ^2 =3C^,s / / — J— 7^-— — , B19 

JO Jo [I + s cos^Oi-\- g cos^02Y'^ 
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y'^dy 



which after obvious substitutions takes the form 

(3cos^6'icos^6'2) = 3C(5f, s) / x^dx i — , 

Jo Jo (1 + sx^ + gy'^Y''^ 

Let s < and g <Q, then 

g = -d? , 
for some non-zero b and d. Which gives for the inner integral in Eq. ( P^ 

y'^dy 



(B20) 



(B21) 



Gi 



(1 — b'^x'^ — d'^y'^y/'^ 
1 1 



y 



d^^yl - b'x' - d'y' d^ 
d 



1 . d?/ 

arcsm- 



Vl — b'^x'^ 



-arcsm- 



dVl - ^'a;^ - d^y^ # v^i - 62a;2 



(B22) 



Therefore 



(3 cos^ di cos^ 62) 



?>C{g,s) \ /•! a;2da; 



ci2 



1 /■! 



Vl — d^ — h^x"^ d Jo 



d 



— / X arcsm 



yr^F 



x^ 



=dx^ . (B23) 



The second integral in the square brackets can be integrated by parts 



Go = - X arcsm ^ 

dJo VI - b^x^ 



dx 



db^ /•! 



x^dx 



1 d 

— arcsm ^ , , , 

d VI - b'^x^ d Jo (l-62a;2)Vl-rf2-62a;2_ 



where 



( 



d 



db^x 



arcsm - 



(B24) 



(B25) 



V Vl-&2x2; (l-62a;2)Vl-t/2_ 523,2 

was taken into account. The last integral in Eq. (|B24|) can be evaluated using the identity 

- h^x^ = -x^ib^x^ ~l)-x^ . (B26) 



Indeed 



b'^x^dx 



1 x'^dx 



x'^dx 



(1 — 62x2)a/1 — d^ — 62^2 Jo \/T^^dF^^Wx^ Jo (1 — 623;2)^X — d^ — 62^2 

(B27) 



Thus 



(3cos^^icos^^2) 



t/2 

1 



1 x'^dx 



1 rf 

arcsin 



VI - t/2 - 62a;2 rf "^ Vr^&^ 

x2dx I 



(l-62a;2)Vl-c/2_ 623,2 



(B28) 
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which gives 



{6 COS fc'i COS 62} = C{g,s)< — 1 — — arcsm- 



62 53 ,/T^^ 

1 d 1 b 1 1 /T. X 

- arcsm^ — -— arcsm^ + — — — > . (B29) 



Finally, 



(3 COS fc'i COS ty2) = C[g,s)< --— — arcsm 



d%^ 63 ./Tzr^ 

which is symmetric in respect to the interchange g <-^ s if both parameters are negative. 
The symmetry breaks if the parameters have opposite signs. Let 

s = b^ , 
g = -d'^ , 

then one has to calculate 

h = C{g, s) 2 / ,^ „ , ,, , : + 



(B31) 



^/T^^WT¥x^ Jo {i + b'^x^)VT^^WT¥x^ 

Id] , , 

- arcsin^ > . (B32) 



The result is 



_ C{g, s) j VI - ci2 + 62 1 _ rf2 
^^ - ~1^ b^ ¥~ ^^ 



6 + VT^rf^+F 



VT^^ 



63 v^r^^ 62C(^,s) d VTT¥j ^ ' 



If s is negative while (/ is positive, the integral is 

C{g, s) \ yi + d^-b^ 1 - 6^ 
^ ~ 62 1 rf2 ^^ 



In 



rf + VI + t/2 - 62 



v^r^^ 



1 ,„<liv^EZ__L__l ,,„„_^\ . (B34) 

rf3 yr^rf2 rf2C((7,s) 6 VTTd^l ^ ' 

For s and g both positive, one can obtain the integral by changing the inverse trigonometric 
functions with inverse hyperbolic functions: 



^(..4-^^mi^.ia..cs,„h- ' 



62rf2 63 VrTrf2 

larcsinh-^ + ^,^^i^^} . (B35) 
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Note that this expression is also symmetric in respect to s *-^ g interchange. The following 
formula valid for all values of s and g, sums up our results 



1 



(cos' e) = — {l + g + s + gs- C{g, s)^l+g + s 



Note, that C{g, s) is given by Eq. (ITO) completely defines the solution for (cos^ 



(B36) 



C. Analytics for perfect coupling 



Analytical solutions for the alignment measure corresponding to perfect coupling of J 
with the axis of major inertia were obtained in Paper I. Here we write down those solutions 
in the form convinient for comparing with the solutions obtained in the main body of the 
present paper. 



For oblate grains, 



a 



3(1+^) 
2s 



1 + s + g 



arcsm 



1 + 9 



2' ^<° 



(CI) 



and 



a 



3(1 + ^) 

4^3/2 

X (ln(l + g) 



-21n(sV2 + (l + s + (7)1/2)) 



2' 



s > 



(C2) 
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Figure Captions 

Fig. 1. Distribution function of J in the grain reference frame for T^s/Ts = 100 

Fig. 2. Zi axis of the external or gas reference frame, XiYiZi, is directed along 
magnetic field. The internal or grain frame, X2Y2Z2, is defined so that Z2 coincides with the 
symmetry axis of the spheroid. 61, ipi, 62, and (/92 are the polar angles in the above frames. 

Fig. 3. Rayleigh reduction factor, R, for oblate grains {g < 0) subjected to a fiux with 
s < 0. A sharp peak corresponds to g = s = —0.5. 

Fig. 4. Rayleigh reduction factor, R, for oblate grains {g < 0) under MHD waves 
(s = —0.5). The alignment is most efficient for flakes {g = —0.5). Due to the intrinsic 
symmetry inherent to the problem the same plot represents the Rayleigh reduction factor 
for flakes {g = —0.5) when s varies from —0.5 to along the x-axis. 

Fig. 5. Rayleigh reduction factor, R, for flakes {g = —0.5) under streaming motions 
corresponding to s > 0. Evidently such an alignment is inefficient. 

Fig. 6. Rayleigh reduction factor, R, for oblate grains [g < 0) when s > 0. The 
alignment measure for prolate grains when s < can be obtained simply by inchanging s 
and g. 

Fig. 7. Rayleigh reduction factor, R, for prolate grains [g > 0) if u_LB (s = —0.5). 
The alignment is marginal. 

Fig. 8. Rayleigh reduction factor, R, for prolate grains [g < 0) and s > 0. This 
situation corresponds to the alignment through streaming along magnetic field lines. 

Fig. 9. The measure of alinment, R, for prolate grains streaming along magnetic fleld 
lines as a function of grain axis ratio y and the velocity ratio w. High degree of alignment 
corresponds to both large y and small w. 

Fig. 10. Rayleigh reduction factor R of prolate grains streaming along magnetic fleld 
lines for y = 10 (upper plot) and for y = 5 (lower plot) for a wide range of the velocities 
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ratio w. It is evident, that for mildly supersonic drift the alignment is marginal. 

Fig. 11. Rayleigh reduction factor R of oblate grains subjected to the drift with u±B 
for J coupled with the axis of major inertia as a result of internal dissipation (upper plot) 
and for negligible internal dissipation (lower plot). 

Fig. 12. Rayleigh reduction factor of oblate grains subjected to drift with u±B. The 
solid line corresponds to the temperature ratio 100, the dotted line corresponds to 10 and 
the dashed line corresponds to 1.1. 

Fig. 13. Rayleigh reduction factor for u±B as a funcion of grain eccentricity for 
prolate grains. The solid line corresponds to the temperature ratio 100, the dotted line 
corresponds to 10 and the dashed line corresponds to 1.1. 

Fig. 14. Rayleigh reduction factor for radiation pressure (s = 100) as a funcion of 
grain eccentricity. The solid line corresponds to the temperature ratio 100, the dotted line 
corresponds to 10 and the dashed line corresponds to 1.1. 



X 




measure\of alignmeht, R 



0.25 



0.2 



0.15 — 



0.1 



0.05 




0.1 



0.1 -0.5 



measure of alignment, R, for s=-0.5 



0.2- 



0.15 — 



0.1- 



0.05 — 



+ 



H \ 1- 



-0.5 -0.4 -0.3 -0.2 

g 



-0.1 



measure of alignment, R, for g=-0.5 



-0.01 — 



-0.02 



-0.03 



-0.04 — 



H \ 1- 



10 



measure of alignment, R 



-0.01 — 



-0.02 



-0.03 — 



-0.04 




-0.5 



measure of alignment, R, for s=-0.5 



-0.01 — 



-0.02 



-0.03 



-0.04 — 



H \ 1- 



g 



10 



measure of alignment, R 



0.22 



0.16 








measure of alignment, R 



0.14 



0.12 — 



0.08 




0.2 



measure of alignment, R, for y=5 and y=10 



0.18 



0.16- 



0.14 — 



0.12- 



0.1- 



0.08- 



0.06 — 




H 1 \ 1- 



0.7 



0.8 






0.6 






^ 0.4 

w 

cd 

CD 



0.2 








-0.5 



-0.4 



-0.3 
s 



-0.2 



-0.1 











"1 I I I I I I I I I I I I I I I I I r 



-0.02 



-1-1 

CD 



•r— I 
1 1 

Cti 

^-l 
O 

^ -0.04 

CO 



-0.06 



J I I I I I I I I I I I I I I I I I L 







8 



10 



g 






•r— I 
1 1 

O 
CD 

P 
W 

cd 

CD 



0.2 







-0.2 



-0.4 








g 



